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Abstract 



It has been argued that the mystery boson X(3872) is a molecule state 
consisting of primarily D^D*^ + D^D*^. In contrast, apparent puzzles and 
potential difficulties have been pointed out for the charmonium assignment 
of X(3872). We examine several aspects of these alternatives by semi- 
quantitative methods since quantitatively accurate results are often hard to 
reach on them. We point out that some of the observed properties of X(3872), 
in particular, the binding and the production rates are incompatible with the 
molecule interpretation. Despite puzzles and obstacles, X(3872) may fit more 
likely to the excited ^Pi charmonium than to the molecule after the mixing 
of cc with DD* + DD* is taken into account. 

PACS number(s): 14.40. Gx, 13.25.Gv. 12.39.Mk, 12.39.Pn 
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I. INTRODUCTION 



The narrow state X(3872) was discovered by the Belle Collaboration [1] in B decay 
and subsequently confirmed by the BaBar Collaboration [2] . Both the CDF [3] and the DO 
Collaboration [4] saw a clear signal of inclusive X(3872) production in pp collision. Since 

X(3872) decays into 'it~^it~ J/ip, it was the most natural to assign X(3872) to one of the 
excited charmonia [5,6]. In fact, the experimental study of the TrnJ/ip mode in B decay 
had been encouraged by theorists, notably Eighteen, Lane and Quigg [7], as a means to 
explore the excited charmonia. However, experiment has revealed unexpected properties of 
X(3872). 

The very narrow width rules out assignment of X(3872) to a natural spin-parity state 
(P = (—1)"^ with C = P) since it would quickly decay into DD. Among the unnatural 
spin-parity states, the 2^Pi(/i^) charmonium of 1"*"" was ruled out by Belle [8] at an early 
stage through the angular distribution of the final 7r"'"7r~ [9]. If we trust the potential- model 
calculations of the charmonium mass spectrum from the past, there is no suitable candidate 
for X(3872) in the close neighborhood of 3872 MeV. Meanwhile, the coincidence of the 
X(3872) mass (3871.9 ± 0.5 ± 0.5 MeV [1]) with the D^D*^ threshold (3871.3 ± 0.5 MeV) 
prompted many theorists [10] to revive the idea of molecular states [11] and to speculate 
that X(3872) may be a loosely bound state of DD + DD* through a color-neutral force. 

Then the Belle Collaboration [12] reported the startling discovery that X(3872) decays 
into ujJ/ip as well, actually 7r"'"7r~7r° off the cu resonance peak, with roughly the same branch- 
ing fraction as 7r~^7r~ J/ip. Since any hadron must have a definite charge parity, the dipion 
TT+TT" in TT^T[~ J /ip ought to be in p-wave {C = (—1)' = —1), namely, the p and its resonance 
tail in / = 1. In the same analysis they saw a signal of the radiative decay mode 'yJ/ip [13], 
reinforcing the C = +1 assignment for X(3872). With the mass difference between D^D*^ 
and D^D* being as large as 8 MeV, isospin violation may be enhanced if a very loosely 
bound molecule-like state should be formed. Then it would nicely explain the large isospin 
mixing in the X(3872) decay. The most recent analysis of the decay angular distributions 
[14,15] favors J^^ = 1++ for X(3872) among the positive C states. Since DD* + DD* 
would be most likely bound in s-wave, if at all, 1"^"*" fits well to the molecule. The molecule 
interpretation has gained steam among some theorists for this reason. If X{3872) is indeed 
a molecule state, it would open a large new field in hadron spectroscopy. Most recently, 
however, a few pieces of crucial experimental information [16], though yet preliminary, have 
appeared against the molecule assignment. We are now in a rapidly moving state of confu- 
sion. 

In this paper we cast more doubt on the molecule interpretation of X(3872). We first 
show that unlike the proton and the neutron inside a deuteron, no long-range potential arises 
from one-pion exchange between D and D . Though it may sound strange, it is a simple 
consequence of the numerical accident, mo* — ^tid — t^i-K — 0, and of the derivative D*D7r 
coupling. This makes the "deuteron-like molecular binding" of DD + DD* highly unlikely. 
Secondly we argue that the observed decay branching of S+ — > K'^X(3872) in B decay and 
the inclusive production cross section of X{3872) in high-energy proton-proton collision are 
both too large for a loosely bound object of binding energy 1 MeV or less. 

We then turn to the charmonium option. We cannot offer a quantitative resolution as 
for the mass spectrum computed in the potential model. It has been acknowledged that the 
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computation involves large uncertainties near and above the open charm thresholds. In the 
1960's the elaborate coupled-channel N/D method was developed in the S-matrix theory 
to compute the spectrum of light hadrons, which were actually the molecule states in the 
present-day language. The N/D method in its simple version iterates the Born amplitudes 
in s-channel. The computation is very close to the coupled-channel potential problem [17]. 
But we were realized at the end how difficult it was to obtain reliable quantitative results 
and could not go much beyond semiquantitative analysis in most cases. The charmonia 
above the open charm threshold share similar uncertainties though the ambiguity coming 
from the high-energy tail is more contained. We leave the difficulty of the X(3872) mass 
with the potential model calculation of charmonia to future study. Instead we focus on the 
production rate and cross section of X(3872). Our estimate shows that because of the very 
loose binding of the molecule, the production rates of the molecule X(3872) should be at 
least an order of magnitude smaller than what we see in experiment. In comparison, the 
charmonium is more easily produced because of a little tighter binding. One obvious obstacle 
to the charmonium option is the large isospin breaking in the decay X{3872) — > p{u!)J/ip. 
However, the u! resonance peak is nearly two full widths outside the phase space boundary. 
We ask if there is any chance that this severe kinematic suppression of the isospin-allowed 
decay into ujJ/ip can bring B{X — > pj B{X — > ujj/ip) up close to 0(1). Our finding 
is that such chance is not ruled out within the small range of experimental uncertainty 
in the mass of X(3S72). Equally or even more serious is the absolute magnitude of the 
X(3872) p{uj)J/ip decay rate; the decay through cc annihilation should not dominate 
over the mode X(3872) — » p{uj)J/il) to hide them. If the charmonium mixes with a fair 
amount of the DD + D*D component, it could generate measurable branching fractions 
into p{Lj)J/ip. The best scenario that emerges is as follows: The X(3872) state is the 1++ 
charmonium that is bound primarily with the gluon-exchange force but mixed to DD +D*D 
of / = through a light quark exchange. This scenario will survive even if X(3872) is above 
the D^D*^ threshold and the decay X(3872) D^D*^ + D°D*^ indeed occurs. We must 
admit that many of the following analysis are only semiquantitative. This is an inevitable 
shortcoming due to the relevant long-distance dynamics that we do not really know in details. 



II. MOLECULE STATE 
A. One-pion exchange potential 

We first show that one-pion exchange produces practically no force between D and D* 
contrary to intuition. In the limit of m/j* = mD+mT^, there exists only a 5-function potential 
in the s-wave channel. Therefore it is incapable of binding D and D . 

Let us define the D*Dt^ couphng by 

Lint = -igD*^rD ■ d^n + h.c. (1) 

The value of g can be fixed by the £)+7r° decay rate to 5'^/47r = 12.8 ± 3.1. A single 

pion can be exchanged only in the i-channel of the DD D*D scattering. (See Fig. 1.) 

For concreteness, let us examine the n'^ exchange potential in D^D*^ D*^D^ . The 
one-pion exchange would be normally the primary source of the binding force. In the rest 
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FIG. 1. The TT-exchange in the scattering DD — > D*D. The intermediate state at the hori- 
zontal dotted Une occurs on mass shell for 7r° exchange. 



of this paper we shall write DD + DD* simply as DD unless it may cause a confusion. 
The potential is extracted from the Born amphtude near the threshold: 



Tb = (e* • q) 



9' 



— — is 



(2) 



where q is the pion momentum, and e and e' are the polarizations of D and D* , respec- 
tively. It is important to notice here a peculiarity of kinematics in the pion propagator. 
Since the decay D* D ir'^ actually occurs with a tiny Q-value, the denominator of the 
propagator is practically zero as compared with in the nonrelativistic limit; 



2 2/ q^ q'^ \^ 2 

-q ~ - (mo* - mo + ) + q , 
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where A = mo* — mo — m^, is very small (~ 7 MeV). Let us denote the static limit of the 
denominator by 



= 2m^A(~ O.lm^). 



(4) 



With the two powers of q from the derivative DD*7c coupling in the numerator, the Born 
amplitude Tg survives only at q^ ^ 2mT^A and varies slowly there. The Fourier transform 
of Tb gives the one-pion exchange potential. Since D*^ T^n^ can occur on mass shell, 
the principal part of the pion denominator is relevant to the potential: 



AmDmD^V{r) = J RerB(q)e*^'--^, 



2 iq r 



/(.■•.q)(e.q)iL^ 



(i^q 



(27r) 



^-i,V.eM(r) + f(e'*.r)(e.r)^^^ 
o 47r r 



+ |^((e*-e)-3(6*-r)(e-r) 



cos jir ^ ji sm /ir 



(5) 
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For the s-wavc potential relevant to binding of DD* into a molecule of l"'"^, only a tiny 
contribution of 0{fi'^) survives aside from the (5-function term since the last term (valid off 
r = 0) in Eq. (5) goes away up to 0{ii'^) after partial- wave projection.^ Contrary to the 
naive expectation, therefore, the one-pion exchange potential is a 5-function in good ap- 
proximation between D{D) and D (D*). Unlike the one-dimensional (5-function potential, 
the three-dimensional (5-function potential cannot generate a bound state when one puts it 
in the Schrodinger equation. The Yukawa potential can arise only from "crossed" multipion 
exchange. Although the dipion-exchange potential of an intermediate range plays an impor- 
tant role in tightly bound nuclei [18] and even in the deuteron by providing a spin-isospin 
independent force, lack of the long range force would not lead to binding of the deuteron. 
Therefore the speculation fails that X(3872) may be an analog of the deuteron on the ba- 
sis of a long range force [19]. Swanson [20] "appended" the one-pion exchange potential 
with a quark-model potential which he calls the "quark Born diagram potential" [21]. He 
turned his quark potential into a meson-meson interaction potential and added to it the 
standard one-pion exchange potential of range I/itItj-, which should be absent between D 
and D according to our argument. For a spatially extended composite state, it is the force 
of the longest range that is most sensitive to binding. It is not clear in his paper how much 
his quark-model potential is responsible for binding of X(3872). While we are unable to 
read quantitative details of calculation in his short paper, it appears that some fundamental 
revision is needed in his calculation of binding, at least in the one-pion exchange part. 

The imaginary part of the Born amplitude Tb due to the on-slicll intermediate state gives 
an absorptive potential: Its s-wave contribution is small: lmV{r) = {g'^n'^ /127rr){e' -e) sin ^r. 
This absorptive potential is nothing other than a tiny channel coupling to ir^D^D^ by pion 
emission from D . The derivative pion coupling or the p-wave emission of the pion in 
D* —>■ Dn is the source of the behavior ImV"(r) ~ //^ as /x ^ 0. 

It is obvious why the one-pion exchange does not lead to the familiar Yukawa potential of 
g-m^r^^. rpj^g pion emission in D* D + tt can occur without violating energy conservation 
so that the emitted pion can travel over long distances. This would otherwise produce a 
long range potential of ~ 1/r and 1/r^, but the derivative coupling erases them leaving only 
the 5-function in the small jj, limit. 

Is it possible to bind D and D through couphng to the coJ/t/j channel? In general, a 
channel coupling can enhance the effective potential for the channel of the strongest force. 
For this channel coupling to occur, D (or D*) must be exchanged. (See Fig. 2.) Since 
the denominator of the D{D*) propagator is ~ 2m£i/£).mtj or larger, the range of this off- 
diagonal Yukawa potential is very short (~ 0.1 fermi). Consequently no strong coupling 
occurs between the DD* and the uJ/ip channel at distances large enough to be relevant to a 
loosely bound molecule. As for the diagonal potential of uJ /ip, the elastic scattering process 
is the so-called "disconnected" quark diagram (Fig. 3) or multi-gluon exchange processes. 
Just as 0(ss) interacts only weakly with tt, p, u) and so forth even at low energies, elastic 



^For the Tr^-exchange in D^D* — D*'^D~, the Yukawa potential appears rather than the oscil- 
latory potential since the on-shell transition does not occur with 77z^*o < m£,- + m^+. However, 
the leading long-range s-wave potential is ~ fPe~^^^/r instead of ~ /x^ cos /xr/r and still totally 
negligible. 
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LoJ scattering is expected to be weak and hardly a source of binding. To summarize, 
we see no chance of generating a bound state in the s-wave DD* channels, no matter how 
the Born diagrams are iterated, by solving the Schrodinger equation or the Bethe-Salpeter 
equation with or without channel coupling, or by approximating the N-function of the matrix 
N/D with the Born amplitudes. 



CO 



D/D* 



D D* 
FIG. 2. The DjD* exchange in DD* ojJ/ip. 




c c q q 

FIG. 3. The disconnected quark diagram oi uJ ^ uJ 

If there is any relevant channel coupling, it would be to the cc channel. The confining 
potential between c and c is stronger than the hadron exchange potentials. The coupling 
between DD and cc occurs through the ■u((i)-quark exchange. (See Fig. 4.) If this channel 
coupling is strong enough, DD and cc can mix substantially. Even when such a mixing 
occurs, the primary source of binding is the confining force of gluon exchange between c 
and c, not the van der Waals force of QCD, i.e. not the hadron exchange force in the DD 
channel. Such a bound state is conceptually not the molecule and dynamically different 
from it . We will later discuss this possibility in more detail. 



B. Decay rate and production cross section 

We compare the observed production rates of X(3872) with the theoretical expectation 
for the molecule state. We are limited to semiquantitative discussion here since accurate cal- 
culation would require much more information of long-distance dynamics and experimental 
input than we have at present. Nonetheless we see a serious difficulty. 
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FIG. 4. The light quark exchange in cc — > DD . 
The X(3872) mass has been given by the four experimental groups as follows [1-4]: 



^X(3872) 



3872.0 ± 0.6 ± O.SMeV (Belle) 
3873.4 ± 1.4MeV (BaBar) 
3871.3 ± 0.7 ±0.4MeV (CDFII) 
3871.8 ± 3.1 ±3.0MeV (DO). 



(6) 



Olsen [8] quotes 3871.9 ± 0.5 MeV as the weighted average from all groups. We should 

compare these numbers with mDO+mD*" = 3871. 3 ±1.0 McV and mD++mD*+ = 3879.4±1.0 
MeV [22]. For rriD and rriD*, we have added the quoted errors ±0.5 MeV of the Review of 
Particle Physics [22] since the errors are correlated between moo and m£,*o and between 
mD+ and mD*+- Although the central values are on the side of mx(3872) > m/jo ± m^j-o, 
even the number from Belle that has the highest accuracy is not conclusive as to whether 
mx(3872) is really above the D'^D threshold or not. We have shown the current ranges of 
"T.x(3872) and niDO + mD*o in Fig. 5. 



3873 



3872 



3871 



3870 
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FIG. 5. The ranges of mx{^872) the value of m£)0 + m£)*o (denoted as L>0 + D*Q). 



In order for X(3872) to be a molecular bound state, X(3872) must be below the D^D 
threshold; mx(3872) > TnDo-\- mD*o. However, Belle apparently saw a large decay branching 
to D^D^* with the preliminary result [16], 



B{X{m2) D°dV)/5(X(3872) ^ t:+t^-J/iP) ^ 10. 



(7) 
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If this result holds, the molecule model ought to be modified into the "virtual state" model 
[23] . For our discussion of the molecule, we proceed with 

?7^x(3872) < ruDO + mc*o = 3871.3 ± l.OMeV. (8) 

This places the binding energy of D^D*^ in X(3872) at less than 1 MeV. 

The production rate of a loosely bound state in decay and scattering is obtained generally 
by convoluting the bound state wavefunction with a production amplitude of its constituents. 
When the range of production dynamics is much shorter than the size of the bound state, 
the production of the constituents and the formation of a bound state factorize. If X(3872) 
is a loosely bound state of D^D*^ , the constituent momenta in the X(3872) rest frame is 



' EBm£)/D*{c::^ 30MeV). The constituents are streaming in parallel with practically no 
relative motion. Even after boosting to the overall cm frame in B decay or in pp colli- 
sion at small rapidity, the relative momentum is still much smaller than the characteristic 
momentum over which the production amplitude A{p,'k) of £)(|p, k)£) (|p — k, s) varies 
significantly. In this case the invariant amplitude of production may be approximated as 

A(p,k)^A(p,0), (9) 

where we have suppressed dependence of A{p, k) on all other variables. Let us express the 
s-wave bound state of DD* in the rest frame with the wavefunction \l/(k) in momentum 
space and the creation operators aj^^ and al,^ of D* and D, respectively: 

/■ / , aLa^ u + con?'. , , d^k , , 

|X(p = 0, s)) = J ^{k) (10) 

For production of X(3872) with momentum p and a given helicity, start with the production 
amplitude of DD with small relative momentum k: 

B/pp ^ D + D* + anything. (11) 

Superpose the production amplitude A{p, k) with the bound-state wavefunction first in the 
cm frame of DD* to get the X production amplitude in the molecule rest frame: 

1 ~ d^k 

=^(k)A(0,k)-— . (12) 



V^EdEd* ' ' ' ' ^(27r)= 

When |k| is so small that ^4(0, k) ~ A(0, 0), we factor out A(0, 0) and integrate over k 
to obtain {Am cm d*)~^^^ "^{0)^(0,0), where ^(0) is the wavefunction at the origin of the 
relative position. Moving back to the overall cm frame, one obtains for the production rate 
of X with momentum p, 

where is the initial four-momentum and the first J2 denotes integration and summation 
over degrees of freedom of all final particles other than X(3872). Because of the very loose 
binding, |^(0)p appears even when D and D* are produced by nonlocal interactions. 
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The wavefunction at origin square | \&(0) p in Eq. (13) gives us a clue about the production 
rates. The value of |^(0)p is fairly insensitive to details of dynamics for a very loosely bound 
state. We estimate as follows: Note first model independently that ^'(k) oc l/(«;^ + k^) near 
the X(3872) pole {k^ = 2EB'mDT(iD* /(^d + ^d*))- The ^(k) of this form is not suitable for 
determining ^(0) since its Fourier transform ^(r) ~ e~^'^ /A-kt is infinite at r = 0. It means 
that the behavior above |k| ~ k matters in determining ^(0). ^(k) must fall off faster at 
large |k| than the simple pole form.^ We regularize ^(k) into the double pole form with 
parameter M, 



where N — is the normalization. Then we obtain 



2 _ kM{M + k) 



|*(0)|^ = ^^— ^, (15) 

For our estimate we choose M ~ 3k for which the regulator term contributes only about 5% 
at r = 1/k in Eq. (14). Then in the ratio to |^(0)^(25)P ~ 0.024GeV^ of ^(25), we obtain 



^(0)x(3872) 



^(0)V>(25) 



Since the wavefunction spreads far out in space, the value at the origin |^(0)x(3872) P is 
quite small. Consequently production of X(3872) is highly suppressed in B decay and in pp 
collision. We proceed with Eq. (16) to compare with experiment. 

B K+X(3872) decay 

For the charmonium production through the decay b — > ccs, the color-suppressed tree- 
interaction is important and the factorization calculation appears to be qualitatively in line 
with experiment. The two-body charmonium decay B~^ — > K~^{cc) has been observed for rjc, 
J/ijj, Xci, and V'(25) [22,24]: 

B{B+ = 9.0 ± 2.7 X 10'^, 

B{B+ ^ K+J/ij) = 10.61 ± 0.15 d= 0.48 x 10'^, 
B{B^ K+Xci) = 5.79 ± 0.26 ± 0.65 x 10"^ 
B{B+ K+'iIj{2S)) = 6.49 ± 0.59 ± 0.97 x 10"^ (17) 

In comparison, only the upper bounds have been obtained for two-body production of Xdo 
and Xc2 that would not occur in the simple factorization limit: 

B(B+ K+Xco) < 4.4 ± 3.3 ± 0.7 x 10"^ 

B(B+ K+Xc2) < 0.09 ± 0.10 ± 0.11 x 10"^ (18) 



^If one proceeds to compute Eq. (12) with the simple pole form for ^'(k), there is a danger to 
overestimate the production rate. 



9 



The Belle Collaboration [25] has determined the product of the branching fractions as 



B{B+ K+X{3872)) x B{X{3872) n+n-J/ij) = 1.31 ± 0.24 ± 0.13 x 10"^ (19) 

The experimental lower bound can be set on B[B^ — >■ i^+X(3872)) with Eq. (19) after 
taking account of the presence of -B(X(3872) Tr~^Tr~Tr'^J/ip): 

B{B+ K+X{3872)) > 2.6 x 10"^ (20) 

However, if the preliminary result of Belle Eq. (7) is used, the lower bound rises to 

B{B+ K+X{3872)) > 1.6 x 10"^ (21) 

Computation of the decay amphtude for S+ — > K'^X{3872) in the case of the molecular 
X(3872) has been attempted, but without a definite numerical result for the branching 

fraction [26]. It is not surprising when one considers uncertainties involved in such compu- 
tations. Nonetheless, if we dare to make a very crude estimate, we would proceed with Eq. 
(13) and compare B{B+ K+X{3872)) with B{B+ K+iIj{2S)): 



B{B+ K+X{3872)) 
B{B+ K+'i/j{2S)) 



A{B+ K+DD* 



A{B+ K+cc) 



2 



X 



^(0)x(3872) 



*(0)V.(25) 



(22) 



where the cc is in 1 and the DD* is in 1++, both near their production thresholds. We 
estimate the first factor in the right-hand side Eq. (22) with the following reasoning. Imagine 
that a cc pair is produced near its threshold, most likely in relative s-wave. When the cc 
invariant mass is below the open charm threshold, the s-wave cc forms a charmonium, J/V', 
ip{2S), r]c, or r]'^. Above the threshold, they form a pair of D^*^ and d''*^ by picking up light 
quarks. Since by assumption the cc pair production amplitude is insensitive to its invariant 
mass, we expect 

A{B+ i^+(L>°D*°)c=+i) ~ 0.5 X A{B+ ^ K+{cc)), (23) 

where the energy Ecc is little below the open charm threshold while Ejy^* is a little above it. 
The factor 0.5 comes from counting of the relevant charmonia and the multiplicity of spin, 
charge, charge parity states of charmed meson pairs after uu or dd are picked up. We use 
Eq. (16) for the second factor in Eq. (22). Combining the two factors together, we reach 
for the molecule X(3872) 

/ fP \3/2 

S(S+ ^ K+X{3872)) ^ 1.3 X IQ'^^-^^^^ B{B+ ^ K+i;{2S)), 

/ T? \3/2 

0.8x10"^ —\ . (24) 

V0.5MeV; ^ ' 

The number in the last line is more than one order of magnitude smaller than the experi- 
mental lower bound ~ 2.6 x 10~^ in Eq.(20). It is more than two orders of magnitude smaller 
if the preliminary result on the D^D*^ mode is taken into account. Physically speaking, it is 
not easy to produce a large composite object like the molecule in the short- distance B decay 
since |^(0)p requires the constituents to come close to each other in the position space. To 
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enhance the production rate to the level of the experimental lower bound, X(3872) must be 
an object of stronger binding. Ironically from this viewpoint, if X(3872) were a bound state 
primarily of D^D*~ instead of D^D*^, the decay branching fraction would be a little closer 
to experiment. 

pp — > X(3872) + anything 

The CDF Collaboration observed 580 ± 100 events of X production in the region of 
rapidity \ri\ < 1 at y/s = 1.96 TeV when X(3872) is identified with 'n~^7r~ J/ip. They are 
about 10% of ipi^S) production events in the same kinematical region [3]: 

ajpp ^ X(3982) + ■ ■ ■)i?(X(3872) ^ tt+tt- J/^) ^ 580 ± 100 

aipp iP{2S) + ■ ■ ■)B{iP{2S) TT+TT-J/iP) 5790 ± 140 ^ ' ^ ^ 

where "• • •" denotes "anything" . The cross section for '^{2S) above is that of the primary 
production ignoring the 10 ^ 20% contribution of the feed-down from the B[B) decay 
[27]. The ratio of the decay branching fractions is 1.6 if X(3872) y^D^D and 0.27 if 
B{X{3S72) D^d\^) = 10 X B{X{3872) tt+tt" J/V^) (cf Eq. (21)). 

a(pp X(3982) + •••) = (0.06 ~ 0.36) x a{pp -0(2-5) + •••)■ (26) 

Wc compare the yields for X(3872) and ip{2S) production using the same argument on the 
amplitude ratio as described for the B decay. 



a{pp ^ X(3872) 



A{pp DD* + 



A(pp — > cc + • • •) 



*(0)x(3872) 



*(0)v-(25) 



2 

X a{pp^'ip{2S) + ■■■). 



(27) 



This leads us to the estimate 



^^^) xa{pp^xl;{2S) + ---) (28) 



The ratio of X(3872) and ip{2S) production is one and half to two orders of magnitude 
smaller than the CDF observation, depending on whether X(3872) indeed decays into D D 
or not. Strong suppression of the DD molecule production should not be a surprise. Such 
suppression was well known experimentally in similar situations: The production cross sec- 
tion of a deuteron [Eb — 2.2 MeV) in pp — > dTr"*" at the cm energy ^/s — 2.98 GeV [28], for 
instance, is only 0.11 ± 0.06 mb, one hundredth of the continuum pp — > priTr"^ production 
cross section, 11.44 ± 0.65 mb, at the same energy. 

To summarize the case of the molecule, lack of the binding force is a serious shortcoming 
on a theoretical side. Despite the numerical uncertainties involved in our theoretical esti- 
mate, experiment shows that production of X(3872) is much stronger than we expect for a 
very loosely bound state. Our conclusion on the molecule applies to the DD* virtual state 
as well since the wavefunction of the virtual state spreads even more than that of the bound 
state. 
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III. CHARMONIUM 



If X(3872) is a charmonium, the most likely candidate is the radially excited ^Pi state 
with 1++. While ^^2(2"+) may look promising, the angular distribution analysis disfavors 
it [14,15]. The charmonium interpretation of X(3872) encounters two immediate problems. 
One is the discrepancy with the potential-model calculation of the mass spectrum [6,5,29] 
and the other is the large decay branching into the 1=1 channel, X(3872) pJ/'^P- We 
do not attempt to propose any resolution for the potential model calculation. Soon after 
the first discovery of charmonia, the mass spectrum below the open charm threshold was 
fitted and reproduced well with the Cornell potential model [32]. Above the open charm 
threshold, calculation must include coupling to the charm-meson-pair channels. While such 
computation was already undertaken even in the early paper of the Cornell model, the 
results involved much larger uncertainties than those below the threshold. A decade before 
the charmonium, strenuous efforts had been made in multi-channel computation of the light 
hadron mass spectrum. However, we were unsuccessful in producing quantitatively accurate 
results and content with the semiquantitative predictions in the strongest-attractive-force 
channel argument.^ With this excuse we will not go into the question as to whether or 
not the channel coupling to DD* and D*D can indeed lower the 2^Pi mass to 3872 MeV 
from the existing predictions of the potential model calculation. As for the production 
in S-decay and pp-coUision, the larger wavefunction overlap of the 1"*"*" charmonium than 
that of the molecule enhances significantly the production rate. However, we arc unable 
to produce quantitatively accurate results, since the charmonium wavefunction is just as 
sensitive to the channel coupling as the mass spectrum is. Furthermore, the approximation 
of the wavefuction-at-origin is less reliable for the charmonium production in B decay. 

Leaving the mass spectrum aside, we instead focus on the issue of the large pJ/ip branch- 
ing fraction. 



A. Kinematical suppression of X(3872) pJ/ip 

Coexistence of the ujJ/ip{I = 0) and pJ/ip{I = 1) decay modes [12] appears as a strong 
argument in favor of the molecule model: 

y--"r-°;y,t'=1.0±0.4±0.3. (29) 

B{X ^TT+TT-J/l/j) ' ^ ^ 

However, this relative branching fraction can be misleading with respect to the magnitude 
of isospin mixing in X(3872). The 7r"'"7r~7r° in the final state 7r~^7r~7r'^J/ip comes from the far 
tail of the lu resonance since ujJ/ip is outside the phase space boundary. In contrast, pJ/ip 



^In retrospect, of course, we were barking at a wrong tree in those days by trying to reproduce 
the quark-antiquark bound states as the "molecular states". Nonetheless, lack of accuracy or 
magnitude of uncertainty in the multichannel calculation of the bound-state spectrum was clearly 
appreciated. The N/D method was the most commonly used technique in those computations. It 
is more sophisticated than the naive nonrelativistic potential model and its extensions [17]. 
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is right on the phase space boundary so that tt'^tt' can come from the lower half of the 
p resonance region. Consequently, tt+tt^tt^ receives much stronger suppression than tt+tt". 
Magnitude of this relative suppression is purely kinematical and sensitive to the X(3872) 
mass even within the small uncertainty of the X{3S72) mass. 

The mass difference nix —irij/^ is smaller than the peak value of the uj resonance (782. 6± 
0.1) by several MeV: 



' -7.5 ± 0.6 ± O.SMeV (Belle) 
-6.1±1.4MeV (BaBar) . 

-8.2 ± 0.7MeV (CDF) 



(30) 



Recall furthermore that mjjo + mD*o > mx(3872) is required to prevent the open charm decay, 
that is, 



"^X(3872) ^ '>TT'D0 + TTljjtO — 3871.3 ± l.OMcV, 



(31) 



which leads to mx — rrij/^ — < —8.3 ±1.0 MeV. The tt+tt resonates about twice the 
half width (ip^ = 4.25 ± 0.05 MeV) below the u peak. Even at the high -mass end of the 
7r"'"7r~7r° invariant mass, the height of the Breit-Wigner resonance is 0.19 ~ 0.23 of its peak 
value for mx(3872) = 3871.3 MeV. Another phase space suppression occurs by the relative 
motion between the "off-shell" cu and J/ip. The s-wave threshold factor |pj/^-,r7r7r| skews the 
Breit-Wigner shape and suppresses the high mass end of 7r^7r~7r°. The combined suppression 
from the two sources is quite severe. (Sec Fig. 6.) According to one computation [31] in the 
molecule model; the observed ratio of Eq. (29) would be reproduced if the effective XJ/i/jV 
coupling ratio \g{XujJ /ip)/ g{XpJ /ip)\'^ ~ is 11.5 ± 5.5. 



physical region 



Phase 




Space 



CO— resonance 



FIG. 6. A schematic drawing of the kinematical suppression near the phase space boundary for 
X(3872) ujj/ij). The p peak sits right on the phase space boundary with the height ~ l/18th 
of the u) peak. In this scale the p resonance profile is almost flat. 

We explore here whether there is a chance to explain the large isospin breaking of Eq. 
(29) by the charmonium decay process, not by the isospin mixing composition of X(3872) 
that is advocated in the molecule model. Assuming that u is entirely responsible for 7r^7r~7r° 
in line with Belle's invariant mass plot, we write out the sensitive part of the decay rate as 



r 

r(X(3872) ^ 7r+7r-7r° J/V^) = Fq / 

^3 



rriuj—Sm 



\V{W)\T^ 



3m. [{W - m^f + F2/4] 



dW, 



(32) 
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where W is the tt+tt^tt^ invariant mass, p(VF) denotes the threc-pion momentum Ptttttt in 
the X(3872) rest frame, and Sm. = m^^ + m,//^. — mx(?,s72) is the distance from the phase 
space boundary to the uj peak. Wc have computed numerically the integral in Eq. (32) in 
ratio to the corresponding quantity for p by varying mx over 3870 MeV to 3872 MeV. In 
carrying this calculation, we have cut off the tt'^tt^tt" invariant mass at 750 MeV, as chosen 
by Belle, and the 7r"^7r~ invariant mass at 450 MeV. The resulting kinematical suppression 
K is put in ratio: 

\ = 13-3 ± 1.7 (mx(3872) = 3871 ±lMeV). (33) 

In order to account for the near equality of the observed pJ/ip and uiJ/ip rates, the production 
amplitude ratio must be such that 



AipJ/ij) 



A{ujJ/ilj) 



= 0.27 ±0.02. (34) 



Being kinematical in origin, this is in general agreement with Reference [31] (0.27^ vs 1/11.3). 
The required number of Eq. (34) is an order of magnitude larger than the electromagnetic 
isospin breaking. However it may not be totally out of line for the isospin breaking due to 
the u-d quark mass difference, which is after all the origin of the 8 MeV mass difference be- 
tween D°(D*°) and D+(D*+). Therefore we next explore the magnitude of isospin breaking 
whether the amplitude ratio of Eq. (34) can be realized by the isospin mass splitting. 



B. Isospin breaking due to mass difference 

The molecule X(3872) state is predominantly made of D^D*^ and mixed with D^D* , 
ujJ/'il) and pJ/ip. The large pJ/ip branching is attributed to the isospin mixed composition 
of DD . In the charmonium X(3872), the particle composition is almost purely in 7 = 
and the isospin breaking occurs during the decay process. Although experiment cannot 
distinguish between them, the two cases are fundamentally different in hadron dynamics. 
They are not different pictures of the same physics related by the quark-hadron duality or 
the like. 

We expect that the main source of the large isospin violation is in the DD* intermediate 
states near the thresholds where the splitting between m^o + m£)*o and + m/)*- has 
the prominent effect. We can make some quantitative discussion of its magnitude with 
the dispersion relation for the X — > p{u!)J/ip decay amplitudes by dispersing the variable 
s = p\. In this way we can at least quantify the magnitude within our limited knowledge. 
Since X(3872) is below the threshold, only the dispersive part exists at s = ^i:(3872)- 

Define the invariant decay amplitude A{s) by 

{V{q)J/^{pr'\{n + m\)X^{m = ^f^f''^''<PyA<l)lP>^Ms) (35) 

for the decay X{P) p/uj{q) + J/ip{p) with P = p + q and s = P^. Wc keep only the s- 
wave term near the DD threshold ignoring the d-wave contribution for our semiquantitative 
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analysis. We write the unsubtracted dispersion relation^ for A{s) integrating along the cut 
which runs on the right-hand: 

TT Jsth S' - S 

The lowest intermediate state is tt+tt^J/V' for pJ/tjj and 7r+7r~7r°J/'0 for (vJ/t/j. The dis- 
persion integral of such intermediate states represents elastic rescattering of p{uj)Jl%l) in the 
final state. In the quark diagrams they are the "disconnected processes" (Fig. 3) that are 
not the most favored processes. The dominant low-energy absorptive parts come from the 
open charm intermediate states DD +DD* which couple more strongly with X(3872) above 
the threshold (s > {mo + mu*)^)- 

Since our primary interest is in the contributions of the D^D*^ and D^D* intermediate 
states of C = -|- to the dispersion integrals, we separate them out as 



TT J {mjjo+mjjtoy^ S — S TT 7(m^++m^«+)2 S — S 

if '-'^'"''••d.-, (37) 

7ri4m2 s' - S 



where the dual sign in front of the second integral is -|- for uJ/il) and — for pJ/il). Define 
s-wave XDD* coupling as approximately isospin invariant as 

Lint = -mxfxDD*DD*''X^, (38) 

and DD* p(oi))J/ijj scattering amplitudes as 

(l^(g)J/^(p)-* I D ipi)D*{p2r') (39) 

X ^e^{qre,{pyeM[iP + Q)xM,{s,t) + (p - q)xM2{s,t)]. (40) 
mx 



The absorptive part for the DD* intermediate state is 

Woo,-+(s) = ^-^^^§^^fxDD*M{s)oo,-+{s) (41) 

with |p(s)| = \J{s — {niD + 'mD*y){s — {uid — mD*Y) /'^rnx and M(s)oo,-+ is the s-wave 
projection of Mi(s, t). The ratio of the / = 1 to the 7 = decay amplitude is given by 

Mi^x)p _ loo - 1-+ ^^2) 



A{rn?x)^ Joo + /_+ + /* 



^We obviously write an unsubtracted dispersion relation since a subtracted one would have no 
predictive power on the quantity of our interest. 
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where 



ds'. (43) 



'{"iD++™D*+)^ s —rrix 
1 /■ ImA(s')oo + ^(s')- 



Since the decay X(3g72) p{uj)J/ip is a long distance process, high intermediate states are 
less important. That is, the absorptive part ImA{s) falls off with increasing s. In order to 
make a numerical estimate, we need to know how far the integrals over s' should be extended. 
Let us choose here the effective cutoff Smax no higher than y/s' — 2mD*{— tud + mD* + ?77.^) 
{i.e., Smax = Sjnax " ('^D + iTiD*)^ — IGcV^ and ~ 0) and see how large the isospin 
breaking can be. In this energy region, the most important s'-dependence is in |p(s')|. We 
approximate the rest of Im74(s') to be constant. In this crude approximation which is almost 
independent of dynamics except for the value of the cutoff of integral, the ratio of Eq. (43) 
takes a simple form particularly when we take the limit oi mo — tud* <^ uid + ttid* and 
mo + mo* — mx < mo* — mo — m^r, 

A{m]^)p _ 7r^mx/2[(mz)- + mD*+ - rrixY^'^ - (m/jo + m£,*o - mxY^^] ^^^^ 
M''^x)uj 2^Sinax " \frn^j2[{m d- + mD'+ - mxY^"^ + {moo + m^.o - mxY^"^] 

For V^max —1 GeV, the right-hand side varies from ~ 0.23 at mx(3872) = ttid + m-jj* to 0.18 
at mx(3872) = '"^D + ~ O.SMeV. These numbers are larger than (m^ — mu)/mj^ ^ 0.05. 
The enhancement arises from the fact that the small number tjid + mD* — mx is made 
less small by the square root threshold factor. Combining this isospin breaking with the 
preceding estimate of the coJ/t/j suppression Eq. (34), we obtain for mx — 3871 MeV (cf. 
Eq. (31)) 

- ^ ^ ^ 2 1.0 ± 0.4 ± 0.3). (45) 

B{X ^Tr+TT-J/ip) \0.2j 13.3 ^ ^ ^ ^ 

We are short of the central value of experiment by factor two with large uncertainties. 
The ratio \A{mx)p/A{mx)u}\ is sensitive to the value of the cutoff Smax- If V^max = 
Y^Smax — {iTT'D + m/)*)^ is lowcrcd from 1 GeV to 700 MeV, the ratio of Eq. (45) would 
become unity. Although choosing the cutoff ^Smax below 2m n* would be unrealistic, taper- 
ing off Im74oo,-+(s') towards ^/s^^ = 2mD* has the same effect. 

We have computed above only the kinematical effect due to the mass splitting between the 
charged and neutral charmed mesons. We must admit that we have stretched the numbers 
to the limit in this estimate. Since we have used the isospin symmetric XDD* coupling, 

however, our (/qq — / h) term does not include the final state rescattering of DD in / = 1 

channel at m'^^^^^2) < < Smax- It can generate additional enhancement or suppression to 
the integrals. But we do not have enough knowledge to analyze such dynamical effects. 

We should learn from the exercise above that although the problem of the large isospin 
breaking is serious, it is too early to reject the charmonium interpretation of X(3872) on 
the basis of the large pJ/ip decay branching alone. 
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C. Magnitude of p{u))J/'tp rates 



We should be equally or even more concerned with the magnitude of the branching frac- 
tions themselves than their ratio. In order to observe the decay mode p{ijj)J /ip in experiment 
at all, it must have a branching fraction large enough to stand out of the annihilation process 
^Pi gggiQQg) hadrons. With the kinematical suppression being so strong, would the 
decays into p{uj)J/ip be still visible ? This problem made some theorists suspicious about 
the charmonium interpretation of X(3872) already at an early stage^ We have no means to 
estimate reliably the magnitude of the coupling Jxdd* nor the amplitude lmA{s)oo-+ in 
Eq. (41). Nonetheless we must address to this question. 

Since we mean by charmonium the cc state that is bound primarily by the confining 
force, mixing of cc with DD is not large in a "clean" charmonium by definition. Therefore 
l^'(O) x(3872)P should not be very far from the value of the simplest potential model that 
ignores the open charm channels. When we take the ratio of the branching fractions in B 
decay to that of the production cross sections in pp collision for X(3S72) and '^{IS), the 
wavefunctions |\E'(0)^(25) P and |^'(0)x(3872)P cancel out. It is a reasonable assumption that 
the ratio of the cc production in s-wave to p-wave in B decay is similar to the same ratio in 
pp collision at low rapidity. If so, we obtain 

B{B K+X{m2)) _ a{pp X(3872) + anything) 
B{B ^ K+^j{2S)) ^ (7(pp ^ ^(25) + anything) ' ^ ' 

Denoting the branching fraction 5(X(3872) tt^tt' J/ip) by b^^+j^-, we have 0.02/6,^+71- for 
the left-hand side with Eqs. (17) and (19) while O.OS/^Tr+Tr- for the right-hand side with 
Eq. (25). The both sides are in fair agreement with each other, which may indicate that 
our crude reasoning on the long-distance physics is not out of line. However, we should be 
concerned with the annihilation decay X(3872) — ggg + qqg. This mode, which is insensitive 
to the D and D* masses, is expected to be two and half orders of magnitude stronger than 
2^ Pi "jJ/ip when the mass is adjusted to 3872 MeV in the estimate by Barnes and Godfrey 
[5]. Since B{X{3S72) -fj/t/j) is about 10% of 5(X(3872) ^ uJ/i/j) according to Belle, 
6^+^- should be quite small: 

b^+^- ~ 0.015 (47) 

with the annihilation mode even in the absence of the D^D*^ mode. For the pure char- 
monium X(3872), Eq. (47) requires with Eq. (19) that ip{2S) and X(3872) arc equally 
copiously produced in 5-dccay in the case of X(3872) y^D^D . X(3872) could be more 
strongly produced than ^p{2S) if X(3872) indeed decays into D^D . This is a serious prob- 
lem: |^'(0)p must be large enough for sufficient production of X(3872) in B decay and pp 
collision while a large value of |^'(0)p potentially makes the annihilation decay too strong 
and the decay X(3872) -k^ ti^ (ji^) J / ip invisible. One escape from this problem is that 
DD is already present in X(3872) and easily decays into X(3872) tt~^tt~{7:^)J/iP by 
the quark rearrangement process. We now look into the DD content of the charmonium, 
namely mixing with DD . 



^See for instance Reference [5]. 
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IV. LARGE MIXING BETWEEN CHARMONIUM AND CHARMED MESONS 



The channel coupling between cc and D^*^D^*^ was studied in the potential model. The s- 
wave channels of cc were studied numerically even in the expanded paper of the Cornell model 
[32] . However, it is not clear how much numerical uncertainty should be attached to the mass 
spectrum involving cc at and above the open charm threshold. Diagonalization must be made 
not simply for the multichannel amplitudes at the energy of a bound state, but for certain 
integrals of them. In the multi-channel N/D method, the tractable approximation close to 
the potential model is to represent the N-function of given J^'^ by the Born amplitudes and 
to set to D{oo) — 1 for the D-function. One advantage of the N/D method over the potential 
model is its notational simplicity in discussing channel coupling. Normalize the partial- wave 
amplitude such that unitarity holds as Imaj(s) = a^j{s)p{s)^^aj{s) and lm[aj(s)^-'^] = — p(s) 
where p{s) is the diagonal phase space matrix [17,33]. Introduce the N and D functions by 
aj(s) = Nj{s)Dj{s)^^ . Suppressing the subscript J hereafter. 



where the rows and columns refer to cc and DD of / = 0. We do not include other channels 
such as ujj/ip and ggg as constituents of X(3872) here since coupling of cc to these channels 
is expected to be much weaker than to DD . We treat the coupling to the weakly coupled 
channels as decay. Such separation would become less clear if their channel coupling were 
stronger. The (2 x 2) D-function is 

OM = /-lr^'''W'''''-'. (49) 

TT Jsth S — S 

The zero of detD(s) gives the mass square of a bound state and the diagonalization matrix 
of D{s) at the zero determines the composition of the bound state. ^ Although it is hard to 
get numerical results in our case, we can make one simple observation about the cc-DD 
mixing. 

We are interested in the possibility that a large off-diagonal element Bi2{s){= B2i{s)) 
for cc ^ DD causes a strong mixing. The Born diagram for cc ^ DD is a light quark 
exchange (Fig. 4). We have already shown in Section 11 that there is practically no force 
between D{D*) and D*{D) in the 1+"^ channel. Therefore B22{s) ~ in the low-energy 
region and D{s) has the pattern of 

where Dij{s) = — J[pN/ {s' — s)]ijds'. It is diagonalized by the orthogonal rotation of angle 
9 that is given by 



^In the case of a single channel Eqs. (48) and (49) combined reduce to solving in effect the cor- 
responding Bethe-Salpeter equation where the kernel is the Born ampUtude. In the nonrelativistic 
limit, therefore, it is equivalent to the potential calculation. 
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where the prime on Dij denotes the first derivative in s. The confining potential of gluon 
exchange for cc <-> cc is no weaker than the quark exchange for cc <-> DD . We can thus set 
a bound on the mixing, 

tan2^<2 ^ |^| < 32°. (52) 

If the mixing of ^ ~ 32° really occurs, the binding force would be enhanced by about 60% 
as a channel coupling effect. 

We make the following observation from this exercise: It is possible for the 1++ char- 
monium to contain a DD component up to one third (~ (tan 32°)^). A stronger mixing is 
possible only if the force of shorter distances in the elastic DD channel should play a role in 
binding. The mass splitting between D^*^^ and should not generate a very large depar- 
ture from isospin symmetry in the DD content of X(3872) since the effective binding force 
is the integral of the Born amplitudes smoothed out over energy, not the Born amplitudes 
themselves at or near mx(3872)- The DD content of X(3872) is approximately in / = 0. It 
should also be pointed out that the binding force due to the channel coupling, determined 
by Di2{m\) is insensitive whether X(3872) is above or below the D^D*^ threshold. 

In the presence of a large mixing, production of X(3872) occurs mainly through the 
dominant cc component. The production is robust since |\[''(0)p is large for the p-wave 
charmonia. On the other hand the DD component is unimportant for production since 
its |^E'(0)p is small. However, the DD* component plays the major role in the decay into 
t:^'k~J/%1) and 7r'^7r~7r'^J/ip since the virtual DD component can decay more easily into 
those channels than cc does. The decay DD p{uj)J/ib is a quark-rearrangement process 
and the strength of the DD binding is unimportant. Since the DD component can make 
up to one third, X(3872) — > p{u})J/'il) is more competitive with X(3872) — > ggg + qqg than 
in the case of the unmixed pure charmonium. In our dispersion relation of the tree-point 
function in Section III, a large DD* component is present when the coupling fxnn* is strong 
and the transition to the DD intermediate state is easy. The 2^Pi charmonium mixed with 
DD was discussed by Meng, Gao and Chao [34] in the case that binding is entirely due to 
the confining force. 



V. CONCLUDING REMARKS 

We have examined the molecule model and the charmonium model for X{3872). The 
main motivations of the molecule idea arc the coincidence of the X(3872) mass with uido + 
m£)*o and the large isospin violation in the decay modes. However, there is no long range 
force to bind D and D into a deuteron-like state. The observed production rates of X(3872) 
in B decay and pp collision are too large for a very loosely bound state. On the other hand 
the charmonium has its share of difficulties; The mass does not agree with the potential 
model prediction of the 2^Pi state and the large decay branching for cc — > ggg + qqg could 
make the experimental signal of 7r+7r^(7r°) J/-?/' hardly visible. One resolution appears to 
be as follows: X(3872) is bound primarily by the confining force between c and c which is 
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boosted by the channel couphng to DD* . Production of X(3872) occurs mostly through its 
cc component. The DD component is in / = in good approximation. The large isospin 
breaking of the decay mode pJ/ip relative to ujJ/ip results from a normal magnitude of isospin 
breaking due to D^*^+ - D^'^> mass difference that is enhanced by the severe kinematical 
suppression of the cuJ/ip mode. This picture is very different conceptually from the molecule 
model or the multiquark model. In this picture, the binding force comes primarily from cc 
and secondarily from the channel coupling. The elastic DD channel provides practically no 
binding force. This is an important distinction from the viewpoint of hadron spectroscopy 
because if our picture is right, X(3872) will not an opening of a flood gate for multiquark 
or molecule states. 

How can we distinguish among the different models and pictures by experiment ? We 
should test the particle content of X(3872). In the case of the charmonium and the charmo- 
nium mixed with DD of / = 0, the production rate is the same for and by isospin 
symmetry as long as it occurs through the dominant interaction b — > cos. Independent of 
dynamics, therefore, we expect for the charmonium X(3872), 

T{B+ K+X{3872)) = r(B° ^ K'^ + X(3872)). (53) 

This equality should hold equally well in the case of the large mixing between cc and DD* 
since X(3872) is produced primarily through its cc component whose |^'(0)p is large. For 
the molecule X(3872) where the DD component has a sizable isospin breaking, the equality 
would be violated since the decay amplitudes into / = and 1, or D'^D and D^D* , are 
dynamical independent. If the molecule X(3872) is made as 

X(3872) = L>°D*° cos a + D+D*~ sin a, (54) 

we can express the decay amplitudes for S+ — > K'^X{3S72) as 

A{B+ K+X{3872)) = Aqq cos a + A+_ sin a. (55) 
Then the decay amphtude for B^ K^X{3872) is obtained by isospin rotation: 

A(B° K°X(3872)) = -A+_ cos a - Aqq sin a. (56) 

The color of the spectator quark (u/d) matches that of the c-quark from b — * cos in B^ — >■ 
K+D^D*° and 5° ^ K°D+D*~, but not in B+ K+D+D*' nor 5° ^ K^D^D*'^. In the 
1/Nc expansion, therefore, A+^ — 0(1/ Nc) x Aqq. Consequently for the molecule X(3872) 
we have 

A{B^ K'^X{3872)) = -(^tana + 0(l/iVe)) x A{B+ K+X{3872)). (57) 

Since D^D*^ dominates over D^D* in the molecule (tan^ a < 1), we predict for the molecule 

r(5° ^ K'^X{3872)) ~ tan^a r(fi+ ^ K^X{3872)) 

< r(5+ ^ K+X{3872)). (58) 

Comparison of the B'^ decay with the B~^ decay clearly distinguishes between the deuteron- 
like bound state of D^D*^ (tan^ a -C 1) and the charmonium. For more general molecules, it 



20 



will determine how much asymmetry exists between the D^D* and the D~^D* components 
inside the molecule X(3872). The BaBar Collaboration recently provided [35] 

B(B° K^X (3872)) /B{B+ K+X(3871)) = 0.50 ± 0.30 ± 0.05, (59) 

which is still inconclusive in distinguishing between the molecule and the charmonium. A 
smaller statistical error will decisively rule out the molecule. 

Although we definitely favor the charmonium over the molecule, even the charmonium 
with mixing still has potential difficulties. In addition to the molecule and the charmonium, 
many other models have been proposed [36]. However, there is less handle to extend our 
semiquantitative discussion to those models. A common difficulty for them is that the pro- 
duction rate is likely much lower than the experimental observation since they arc generally 
objects of large spatial spread. The branching fraction ratio of Eq. (59) is useful in distin- 
guishing among them. This ratio is unity for the hybrid (ccg) and the glueball {gg). It is 
far away from unity for the cucu. 
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